Quantum Fluctuations and Dissipation in Thin Superconducting Wires 
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We investigate quantum fluctuations in thin superconducting wires. We demonstrate that quan- 
tum phase slips dominate the system behavior at low temperatures and are well in the measur- 
able range for sufficiently thin wires. We discuss the effect of dissipation, predict a new quantum 
superconductor-to-metal (insulator) phase transitions for wires with thicknesses in the 10-nm range, 
evaluate the resistance R{T) of such wires and compare our results with recent experimental findings. 



It is well known that fluctuations wash out the long- 
range order in low dimensional superconductors . Does 
this result mean that the resistance of such superconduc- 
tors always remains finite (or even infinite), or can it 
drop to zero under certain conditions? A lot is known 
about the behavior of two-dimensional (2D) supercon- 
ducting films where physics is essentially determined by 
the Kosterlitz-Thouless-Berezinskii (KTB) phase transi- 
tion 1^. In quasi-lD superconducting wires below the 
(mean field) critical temperature Tc a nonzero resistivity 
can be caused by thermally activated phase slips (TAPS) 
1^ . This effect is of practical importance at temperatures 
close to Tc where the theoretical predictions have been 
verified experimentally |^. However, as the temperature 
is lowered the number of TAPS decreases exponentially 
and no measurable resistance is predicted by the the- 
ory 1^ at T not very close to Tc- Nevertheless, the ex- 
periments by Giordano clearly demonstrate a notable 
resistivity of ultra-thin superconducting wires far below 
Tc- More recently strong deviations from the TAPS pre- 
diction in thin (quasi-) ID wires have been also demon- 
strated in other experiments 

The natural explanation of these observations is in 
terms of quantum fluctuations which generate quantum 
phase slips (QPS) in ID superconducting wires. How- 
ever, first estimates for the QPS tunneling rate derived 
from the time-dependent Ginzburg-Landau based theo- 
ries turned out to be by far too small to explain the 
experimental findings jsj (see ||^ for more details). 

More recently, the present authors |9| developed a mi- 
croscopic theory describing the QPS phenomenon and 
demonstrated that in sufficiently thin wires QPS effects 
are well in the measurable range and may lead to a 
nonzero wire resistivity even at T = 0. Moreover, the 
existence of a new superconductor-to-metal (insulator) 
phase transition as a function of the wire thickness was 
pointed out in 

In a present paper we extend our theory Q in sev- 
eral important aspects, in particular providing a more 
detailed discussion of the QPS action in various limits 
and paying attention on dissipative effects outside the 
QPS core. We also discuss a possible explanation of the 
recently observed negative magnetoresistance [10[ within 



the framework of our QPS scenario. 

The model. Our calculation is based on the effective 
action approach for a ECS superconductor The 
starting point is the partition function Z expressed as 
an imaginary time path-integral over the electronic fields 
-0 and the gauge fields A, with Euclidean action 
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Here (3 ^ 1/T, $(V) = -V^/2m - fi, en.i denotes the 
background charge density of the ions, and h — ks = 1- 
A Hubbard-Stratonovich transformation introduces the 
energy gap A as an order parameter and the electronic 
degrees of freedom can be integrated out. What remains 
is an expression for the partition function in terms of 
an effective action for A, V and A, with a saddle-point 
solution lAI = Ao and V = A = 0- We obtain 
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where the superfiuid velocity = {l/2m)[Wip — 2eA/c], 
the chemical potential for Cooper pairs ^ — V — (p/2e, 
and A = lAle*"^ have been introduced. 

Effective action for QPS- The effective theory is con- 
structed by expanding up to second order around the 
saddle point in $ and Vs to obtain the electronic po- 
larization terms A phase-slip event in imaginary 
time involves a suppression of the order parameter in 
the phase slip core (i.e. inside the space-time domain 
^ ^ 2;o, T < To), and a winding of the superconducting 
phase around this core. The total QPS action Sqps can 
be presented as a sum of a core part S'corc around the 
phase slip center for which the condensation energy and 
dissipation by normal currents are important, and a hy- 
drodynamic part outside the core 5'out which depends on 
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the hydrodynamics of the electromagnetic fields and dis- 
sipation due to the presence of quasiparticles above the 
superconducting gap. 

In what follows we will consider sufficiently thin wires 
with the cross section S" < A|, where Al is the Lon- 
don penetration length of a bulk superconductor. Due 
to scattering on impurities and boundary imperfections 
the electron mean free path I in such wires is typically 
much shorter than the coherence length of a clean sample 
I <C ^0 = vf/2A. Here we restrict our attention to this 
physically important diffusive limit. Assuming that out- 
side the QPS core the magnitude of the order parameter 
field is not suppressed |A| = Aq for S'out we obtain [p|Jl2| 
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where the integration runs over |x| > xqJtI > tq. In 
general the kinetic inductance L and the kinetic capac- 
itance C in (|l]) depend on frequency u and the wave 
vector k | |T^ . In the limit of low uj and small k we 
have L = AirXj^/S and C ~ Se'^NoUs/n, where Hs and n 
are respectively the superconducting and the total elec- 
tron density. In (|l|) we also introduced the capacitance 
C = Se'^NoTin/n which we will drop from now on in the 
limit Us ^ Hn = n ^ Us at low T. 

The geometry and screening by dielectrics outside the 
wire are accounted for by the capacitance per length 
C and the inductance times length L that replace the 
g2 _|_ B2_j-gj.j^_ jpqj. wires transverse screening is ir- 
relevant and we can retain only one component of the 
vector potential]. The expressions for C and L also de- 
pend on the relevant space and time scales as well as 
on the wire geometry. In the ideal case of a cylindri- 
cal uniform wire for fcro <C 1 (tq is the wire radius) one 
has C = er[21n(/fcro)]~-^ and L — 2 ln(/fcro)/c^, c is the 
speed of light and the dielectric constant of the sub- 
strate. In practice the details of the wire geometry can 
be very complicated (e.g. the cross section S is not con- 
stant along the wire, i.e. the wire is never uniform) and, 
on top of that, other (metallic) objects can be located in 
the vicinity of the wire. The above effects lead to an effec- 
tive cutoff of the logarithmic dependence on k at the scale 
k ^ 1/d with d depending on experimental details (e.g. 
d can be a typical scale of the wire inhomogeneity or the 
distance to the metallic groundplane). Here we will stick 
to a simplified model and assume C = er[2 ln((i/ro)]~^ to 
be constant at all relevant distances. As to L, its particu- 
lar form turns out to be unimportant for thin wires with 
\/S < Xl in which case the kinetic inductance always 
dominates L ^ L. In addition to the above kinetic and 
electromagnetic effects the expression (|l|) accounts for 
dissipative currents outside the core. The corresponding 
contribution is described by the last term in Eq. (|^). 



As to the core contribution, it consists of two terms 
5.o.e = ^iVoAg^roXo + I J2 ^\Eic.^^)f. (2) 

The first part is the condensation energy that is lost in- 
side the core and the second part defines the energy of 
dissipative currents in the core during a phase slip event. 
Here a is the normal state conductance of the wire: we 
already made use of the fact that the typical QPS fre- 
quency is sufficiently high ||^ I/tq ^ Ao, therefore dis- 
sipative currents inside the core are insensitive to super- 
conductivity. It is also important to emphasize that no 
gradient terms for A (both in space and in time) should 
be added to Such terms can be recovered only by 
expanding the effective action in powers of to and k. For 
fast processes (like QPS) this expansion becomes obvi- 
ously incorrect and it is necessary to carry out a more 
careful treatment of polarization terms in the action. For 
the QPS event with I/tq > Aq this treatment yields |l^] 
b - ln(I/(2AoTo) + C'^/xl), where ^ = ^JD/21\q is the 
coherence length of a dirty superconductor, D — vpl/S. 

Variational procedure. In order to evaluate the QPS 
action Sqps ~ Score + S'out wc wiU usc a variational ap- 
proach which consists of several steps. 

We first minimize the hydrodynamic contribution S'out 
with respect to the potentials V and A. As a result we 
arrive at the saddle point conditions which link the po- 
tentials to the phase variable outside the core. Making 
use of the fact that for thin wires one has L ^ L, C ^ C, 
we obtain in the Fourier representation 
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With the aid of (|,|) one can rewrite the action Sout only 
in terms of the phase variable ip{T,x). Then minimizing 
this part of the action with respect to ip and keeping 
in mind the identity dxdr^p — drdx^ = 2ttS{t,x) (which 
follows from the fact that after a wind around the QPS 
center the phase should change by 2n) we find 
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The general expression for the function Q in (|^) is some- 
what tedious and is not presented here. In the following 
limits substantial simplifications can be achieved: 
i) w ^ a{u!, k)Sk'^/C. The function (||) has the form 
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ii) a{Lu, k)Sk^/C < w < a{uj, k)Sk^/C. We find 
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iii) u! <^ cr(w, k)Sk'^/C. The function Q again acquires 
the form (||) but with C substituted by C. 

As a last step of our variational procedure we minimize 
the total QPS action 5coio + S'out with respect to the core 
parameters tq and xq. Let us first neglect dissipation by 
formally putting ct = in (^) and (||). Then solving the 
equations 

dSgps/dro = 0, dSqps/dxQ = 

and treating & as a constant (i.e. neglecting its weak 
dependence on tq and xq) we obtain 
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where cq = 1/ V LC is the velocity of the Mooij-Schon 
mode [ p^ which determines the space-time asymmetry 
of the core. The total action for a single QPS reads 
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where = X'^ + c^(3^, X is the wire length and 
/i = (7r/4e2)y^ C/L. The first term in (||) represents 
the core action S'corc, the second term defines S'out (for 
simplicity we chose the cutoff by integrating outside 
the ellipse {x / xq)"^ + {t / tq)"^ > 1). Substituting L — 
47rA|/5 = l/2ne^N„AoDS into (|[) at T < Aq we find 
xo — Co To — n^/^/b, i.e. the core size xq is of the order 
of the superconducting coherence length f , and the QPS 
time To ~ ^/co <^ I/Aq. This result justifies the above 
conjecture that the typical QPS frequency is higher than 
Ao and demonstrates why our core action S'corc oc xqTo is 
much smaller than that found within the TDGL analysis 
1^ which yields the QPS frequency of order Aq. 

Let us now include dissipation. At high frequencies 
dissipative currents flowing both inside and outside the 
core are important and should be taken into account even 
at T = 0. The dissipative contribution from S'out is ob- 
tained from (^), (^. After a simple integration one finds 



aS/e^xo- 



(10) 



This expression is nothing but the Caldeira-Leggett dis- 
sipative action of a normal conductor with the cross sec- 
tion S and the length ~ .To. Similar expression defines 
the dissipative contribution from the core S^^f.^. 

If a is small one can treat dissipative terms perturba- 
tively. This is sufficient as long as ITS' < e^/i^. It is easy 
to check that in the practically important Drude limit 
a = 2e^NQD the above condition would mean ^ > cq/Aq. 
This condition is never satisfied for realistic parameters. 



Therefore in this limit dissipation cannot be treated per- 
turbatively and our variational procedure should be mod- 
ified. Under certain simplifying assumptions one can find 
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where r — cqTo/xo- The strong dependence of (|TT|) on r 
enforces the minimum condition r « 1, i.e. the asymme- 
try of the core remains approximately the same as in the 
underdamped limit. Under this condition the whole ac- 
tion Sgpg + Sqpg can be easily minimized with respect 
to Xq and we obtain [0 « cqTq ~ £,^/a and 
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Metal- Superconductor phase transition. The next step 
is to consider a gas of QPS's in a superconducting wire. 
We also assume that an applied current / (much smaller 
than the depairing current) is flowing through the wire. 
Substituting the saddle point solution ip — (p(x — 
Xi,T — Ti) into the action and keeping track of the addi- 
tional term ^ dr J dx{I /2e)dx'p |l^, we find 



Sn = nafj, — fi UiVj In 



( PlL 



^/^Z^.T, . (13) 



The quantity pij = (co(Ti — Tj)^ + [xi — Xj)^)^/^ defines 
the distance between the i-th and j-th QPS in the (x, t) 
plane, vi — +1 (—1) are the QPS (anti-QPS) "charges", 
and $0 = hc/2e is the flux quantum. Only neutral QPS 
configurations with vtot = X^T — ^ (and hence n even) 
contribute to the partition function 

For / = Eq. ( |l3| ) defines the standard model of a 2D 
gas of logarithmically interacting charges Vi. The effec- 
tive (small) fugacity y of these charges is 



y = xoToSexp(-a^) 



(14) 



where B is the usual fluctuation determinant which we 
roughly estimate as -B ~ a^/xoTQ. From the Coulomb 
gas analogy, we conclude that a KTB phase transi- 
tion 1^ for QPS's occurs in a superconducting wire at 
/i = /i* = 2 + 4TTy w 2: for /i < /i* the density of free QPS 
in the wire (and therefore its resistance) always remains 
flnite, whereas for > ^^* QPS's and anti-QPS's (AQPS) 
are bound in pairs and the linear resistance of a super- 
conducting wire is strongly suppressed and T-dependent. 
We arrive at an important conclusion: at T — a ID 
superconducting wire has a vanishing linear resistance, 
provided the electromagnetic interaction between phase 
slips is sufficiently strong, i.e. /i > /i*. 

The above analysis is valid for sufficiently long wires. 
For typical experimental parameters, however, X < cq(3 
(or even X <^ cqP), and the finite wire size needs to 
be accounted for. Here we consider the physical situa- 
tion with nonvanishing normal conductivity of the wire 
cr = cTqp (and thus nonzero dissipation) even far from 
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the the QPS core. This situation can be reahzed in the 
presence of quasiparticles above the gap due to finite tem- 
perature and/or nonequihbrium effects. In this case our 
consideration should be modified as follows. 

We first apply the 2D scaling |^ diy — {2 — and 
— — 47r^/i^?/^, where ^ and y depend on the scaling 
parameter I. Solving these equations up to / = Zx = 
\n{X/xo) we obtain the renormalized fugacity y — y{lx)- 

For larger scales I > Ix only the time coordinate mat- 
ters. At sufficiently low frequencies the inter-QPS inter- 
action is determined by the function (|7|) and the problem 
reduces to a that of a ID Coulomb gas with logarith- 
mic interaction. Therefore, (for y ^ 1) further scaling is 
defined by diy = (1 — ^)y and 9;7 = 0, where 

7 — TTSaqp/2e^ is the dimensionless "quasiparticle" 
conductance of the wire. For 7 > 1 the fugacity scales 
down to zero, which again corresponds to a supercon- 
ducting phase, whereas for 7 < 1 it increases indicating a 
resistive phase in complete analogy to a single Josephson 
junction with ohmic dissipation. The phase transition 
point again depends on 5", but also on the wire length X 
and the value Uqp (see below). 

Wire resistance at low T. At any nonzero T the wire 
has a nonzero resistance R{T, I) even in the "ordered" 
phase fi > fi* (or 7 > 1). In order to evaluate R{T) 
in this phase for a long wire we proceed perturbatively 
and first calculate the free energy correction SF due to 
one bound QPS-AQPS pair. The one QPS-AQPS pair 
contribution SF to the free energy of the wire is 



of the free energy reads 
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where p = (cqT^ -|- x'^Y^^. For nonzero / the expres- 
sion in Eq. ( [l5|) is formally divergent for (3 ^ 00 and 
(after a proper analytic continuation) acquires an imagi- 
nary part Im SF This indicates a QPS-induced instabil- 
ity of the superconducting state of the wire. The corre- 
sponding decay rate F = 2lmSF defines the total voltage 
drop V across the wire (see for more details) . For the 
wire resistance R{T, I) — V/I this yields R cx T^t^-^ and 
R cx /^A'-s foj. y ^ (j,^/ and T < $0-^ respectively. For 
thick wires with fi > fi* , we expect a strong temperature 
dependence of the resistivity. For thinner wires the tem- 
perature dependence of the resistivity becomes linear at 
the transition to the disordered phase in which our anal- 
ysis is not valid. At T ^ ^al/c we expect a strongly 
nonlinear I-V characteristic y ~ in thick wires, and 
a universal J^(/i*) = 2 in thin wires at the transition into 
the resistive state with V ^ I, i.e. f = 1. Note that 
in contrast to the KTB transition in 2D superconducting 
films, the jump is from v = 2 to 1, instead of = 3 to 1. 

For a short wire X < co/T we again proceed in two 
steps. A 2D scaling analysis yields the "global" fugac- 
ity y. In analogy with the resistively shunted Josephson 
junction uM, the voltage drop from the imaginary part 
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giving R (X 2^27-2 ^ ^ 727-2 respectively at high 
and low T. Here tq is defined from the high frequency 
cutoff in (^: fo ~ XC/e'^j. The above result is valid for 
7 > 1 and also for smaller 7 at not very small T . At 
T ^ in the metallic phase the resistance becomes |ll[ 



R = Saqp/X 



(16) 



i.e. R is just equal to the quasiparticle resistance of the 
wire whereas the superconducting channel is blocked due 
to quantum fluctuations. 

Discussion. Let us compare our predictions with ex- 
perimental results P,|^,[l0|] . Taking - 10 nm and 
Cr = 1, for typical system parameters kp^ ^ 0.2 nm 
< / - 1 ^ 10 nm < ^ ~ 10 nm < ^0 Al - 100 
nm we obtain the velocity cq/c = cms/c ~ (v^/lOAi), 
/i ~ 30(\/S'/Al) and a 5-^ 10. This estimate yields the 
core action 5'core — ap, \0 in agreement with 

For the quoted parameters, we predict the supercon- 
ductor to metal transition at a wire thickness ^fS ~ 
Al/15 < 10 nm. This prediction also agrees with the 
results of Giordano, who finds that wires with tq — 
\fSpK K, 8 nm have a resistivity that saturates at a mea- 
surable level at low T, whereas the resistivity of thicker 
wires |^] tq > 13 nm always decreases with T. 

Another remarkable feature is that the classical-to- 
quantuni crossover temperature T* was found to be quite 
close to Tq for sufficiently thin wires ||]. Compar- 
ing our quantum action 2S'core with the classical expo- 
nent Q we immediately arrive at a simple estimate for 
T* K, I^q"^ Cq"^ I ^^"^ . For the above parameters it yields 
T* ~ 10A(T*), i.e. for thin wires one indeed expects this 
crossover to happen quite close to Tc- 

Independent measurements of i?(r) for superconduct- 
ing wires have been carried out in Refs. ||^,p^, where sys- 
tematic deviations from classical predictions |^] for thin 
wires have been also reported. Although the overall trend 
|^,|o| is similar to that observed in Q the shape of some 
experimental curves look quite different from It was 
argued that these quantitative differences are due to 
granularity of the wires used in the experiments |^ . How- 
ever, the variations of S were reported to be moderate 
in 1^. If so, this experimental feature can only cause 
a somewhat nonuniform distribution of QPSs along the 
wire because the QPS fugacity increases with decreasing 
S . With trivial modifications our theory can be applied 
to this situation as well. 

In order to proceed with our comparison let us recall 
that the wires in were quite short X ^ \ ^ 2/im 

whereas the wires investigated in Q were typically by 
one or even two orders of magnitude longer. At not very 
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low T one can put aqp ~ a and estimate the parame- 
ter to be of order 1 K or even bigger for the sam- 
ples p|,[Tot. At the same time for the samples ||[ is 
typically below 10 mK. Thus the difference between the 
results 1^ and [pHlOl] can be attributed to the different 
behavior of the function Q (||) for different frequencies: 
the form (||) should be applied in the case [|j whereas 
the samples in the interesting temperature range 

should be rather described by means of (|^). E.g. it ap- 
pears that the data of fig. la of jl^ can be (at least 
qualitatively) described by the dependence 

R{T) cc fiTfof-'-^ (17) 

both for thicker 7 > 1 and thinner 7 < 1 wires. The re- 
sistance of the latter - according to our analysis - should 
increase with decreasing T. This is just what has been 
found in |l^. The crossover from this behavior to that 
with decreasing R{T) for thicker wires can be interpreted 
as an indication of the phase transition at 7 = 1. 

Another interesting feature to be discussed is the neg- 
ative magnetoresistance of the wires observed in [|l^ . At 
the first sight this feature contradicts to our QPS sce- 
nario: in a (sufficiently strong) magnetic field H the gap 
Aq is partially suppressed and the barrier for QPS should 
decrease. Hence, the fugacity y and the wire resistance 
i? - in contrast to jl^] - should increase with H. How- 
ever, if one includes dissipative effects outside the core 
into account the picture can change drastically. Indeed 
aqp — aun/n strongly depends on the relation between 
T and Aq (T, H) . At sufficiently low T a decrease of A 
may lead to an exponential (aqp oc exp(— Aq/T)) increase 
of the number of quasiparticles and - therefore - dissi- 
pation. Thus we have two effects: the field H increases 
the QPS fugacity y but also increases dissipation 7 which 
suppresses quantum fluctuations. It is quite obvious (e.g. 
from Eq. (|l^)) that the second effect may dominate in a 
certain parameter region and the resistance will decrease 
with increasing H. At very large H the gap Aq will be 
suppressed and the resistance will increase back. This 
reentrant behavior was observed in . 

Note added. A clear experimental evidence for a dis- 



sipative phase transition controlled by the normal state 
wire resistance R was obtained very recently in ref. [ p^ . 
In these experiments several uniform wires with thick- 
nesses in the 10 nm range were studied. The wires with 
R < Rq — 7r/2e^ ~ 6.5 Kfl were found to go supercon- 
ducting at low T while the resistance of the wires with 
R > Rq showed no tendency to decrease with tempera- 
ture even well below Tc- It appears, therefore, that the 
results |l^ can be interpreted in terms of the dissipative 
phase transition at 7 = Rq/R = 1 predicted in our pa- 
per. We are grateful to A. Bezryadin and M. Tinkham for 
communicating their results to us prior to publication. 
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